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Abstract 

Liquid-gas phase transition in statistical mechanics is a long-standing 
dilemma not yet well explained. In this paper we propose a novel ap- 
proach to this dilemma, by: 1). Putting forth a new space homogeneity 
assumption. 2) . Giving a new formulation- the mean distance expansion, 
instead of Mayer's Cluster Expansion, to calculate the intermolecular po- 
tential and partition function for a classical system. 3). Explaining how 
the separation of two phases occurs below the critical temperature T c and 
what is the gap between two phases. 4). Calculating the physical quanti- 
ties in a system of coexistent vapor and liquid, and comparing them with 
the experimental results. Qualitative, and some quantitative, consisten- 
cies are obtained. So the statistical explanation on first order liquid-gas 
phase transition is solved in principle. 

PACS numbers: 05.20. Jj, 64.60.-i, 64.70.Fx 

Gas-liquid phase transition is a fundamental problem remaining unsolved for 
centuries in statistical mechanics since the VDW equation of state was developed 
in 1873. Since then, a great deal of attention had been focused on providing 
statistical explanations for the separation of phases by many physicists [1-1 7]. 
This letter provides a new way to solve it. Based on space homogeneity as- 
sumption and mean distance expansion, instead of Mayer's Cluster Expansion, 
to calculate the intermolecular potential and partition function for a classical 
system, this work at first time reveals the existence of the energy gap between 
two phases; offers a rational physical explanation on phase transition and latent 
heat; and demonstrates how to calculate physical quantities in a system of co- 
existent vapor and liquid: such as densities of coexistent gas and liquid, vapor 
pressure and latent heat, etc. As examples, we discussed gas-liquid transition 
for both the VDW fluid and the modified VDW fluid. The experimental results 
provide qualitative, and some quantitative supports to the new theory. 

H.Eyring and D.Henderson etc, mentioned their space homogeneity assump- 
tion (it is referred as EHSE SHA in this paper) in their book[18]: "we shall make 
the unfounded assumption that the molecules are randomly distributed through- 
out the volume and the number of molecules lying near a particular molecule in a 
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concentric shell ranging from r to r + dr is dN = N/V4irr 2 dr, with r > cr. "They 
derived the VDW equation of state from their EHSE SHA, but the total poten- 
tial energy and the partition function are independent of the mean distance of 
molecules in their theory. 

Now we put forth a new Space Homogeneity Assumption (SHA): In a clas- 
sical system in equilibrium with volume V, particle number N, when gravity 
and boundary effects are omitted, each particle averagely occupies volume V/N. 
We will demonstrate in detail later that the SHA will exhibit the correlation 
between the mean distance of molecules and the partition function. 

Both the SHA here and the EHSE SHA should be regarded as independent 
theoretical assumptions, which cannot be proved by statistical mechanics. Their 
justifications lie in the agreement between results derived from them and the 
experimental facts. Now we reevaluate the partition function under the SHA. 

We assume that the potential energy u{rij) between particle i and particle 
j is shared by these two particles together and is divided between both equally. 

JV-l 

We call Ui — \ u { r ij) the potential belonged to particle i if only the pair 

interaction exits. According to the SHA, the volume occupied by particle i 
is equal to V/N = 1/n = r 3 , where n is the particle number density and r is 
the mean distance of molecules. For convenience, we presume that this volume 
is a sphere with diameter d = {(S/irn) 1 / 3 = or in a three-dimensional space, 
and particle i is at its center. We then draw many spherical shells centered 
by particle i outside this sphere. The thickness of each spherical shell is d. 
From inside to outside, we number them S\,S2 and so on in turn. Then the 
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volume of the shell Sk is Vk = f 
the particles in this shell is Nk — nVk- The potential energy of particle i is 
Ui = 1/2 ^2 Nku(rk), where u{rk) is the potential contribution of a particle in 

k 

the shell S k to particle i, (k - l/2)d < r k < {k+l/2)d and Y, N k = N - 1. The 

k 

total potential of the system is U = NUi. The potential energy u{rk) can be 
expanded at the location rk = kd, which is the distance between the midpoint 
of the shell Sk and particle i: u{rk) = u{kd) + u'{kd)(rk — kd) + • ■ •. If only 
taking the first term of this expansion, i.e. u{rk) ~ u(kd), we have the partition 
function 
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with A thermal wavelength. The Eq.(l) clearly indicates that the potential (as 
does the partition function) of the system is correlative with the mean distance. 
Such a relation plays a vital role in the discussion of the phenomenon of liquid- 
gas coexistence. When in Eq.(l) taking the Sutherland potential: u(r) = oo if 
r < a, and u{r) = — s{a/r) 3 if r > a, where a is the diameter of the hard core 
and e is a constant, we obtain the partition function of VDW fluid: 
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with a = N k e^, b = |ttct. The VDW equation is P = k B T-J^ (lnZ ) 

k 

The chemical potential of VDW fluid can be obtained by 
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Letting 
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, , 3, Wirma 
/i= ^T +ln6+ 2 ln "W 



and substituting Eq. (3) into (4), we have 
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with reduced quantities T* = T/T c and n* = n/n c . The Eq.(5) is a function 
of reduced quantities only. Now we can define the reduced chemical potential 
~p* = ~p/~p c . When T < T c , there are three numerical values of n corresponding 
to a given chemical potential fi, which forecasts the coexistence of multiple 
phases. 

For a coexistent liquid-gas system we have the conditions of equilibrium: 
P (T, n g ) — P (T, ni) and \x (T, n g ) — /i (T, m). Thus the solutions of these two 
equations, n g = n g (T) and n; = ni (T) (the densities of the coexistent vapor 
and liquid, respectively), can be obtained by a numerical method, as shown 
(dashed line) in Figure 1A. This curve is equivalent to that given by Maxwell 
rule in thermodynamics. 

The dashed curve for VDW fluid in Fig.lA can be described by the empirical 
formulae (firstly introduced by Wei Wang ) 
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Formula (8) suggests the critical exponent [3 = 0.5 when T* -» 1. Formula 
(9) states the well-known law of the rectilinear diameter. Comparing with the 
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Figure 1: Comparison between the experiment and the theoretical results in 
this paper: A) Densities of coexisting vapor and liquid. B) Vapor pressure. C) 
Latent heat, solid line-experimental curve from Rcf. [20-22], dashed line — VDW 
fluid, dot-line — modified VDW. 



numerical values of the solutions , the inaccuracy of the formula (6) is less than 
0.1% when T* > 0.6, and the inaccuracy of the formula (7) is less than 0.1% 
when T* > 0.71. The percentage of inaccuracy increases as the temperature 
drops below T* = 0.71. 

Substituting them into VDW equation, we then get the vapor pressure P = 
P(T). The reduced vapor pressure is expressed as 

p * = ^T*_ _ 2 = 8r£T _ 

3-n* 9 3-n* 1 ' 

as shown in Figure IB (dashed line). Substituting /i* into Eq. (10), we obtain 
another expression of vapor pressure 

p * = 3nf exp [jW (n* - n* g ) - ^ (nf - nf)} - nf 
l-exp^K-n^-g^K 2 -nf)] 

This relation can be expressed by the empirical formula 

In P* = 3.5204 - 3.5660/T* (12) 

There is a gap of internal energies between the two phases. From Eq.(2) 
the potential energies per particle in vapor and liquid are Ui g — —an g = 

-l/2J2^ku(kar^) and Uu = -ani = -\/2^N k u{koirl) respectively. The 

fe fe 
difference between the two is dU g i = Ui g — Uu = a(n; — n g ), which is the gap 
of internal energies between the two phases due to the different mean distances, 
as shown in Fig. 2, where U* = U/U c . The average distance of gas molecules 
is larger than that of liquid molecules. The interactions between gas molecules 
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Figure 2: The relation between the potential energy per particle and densities. 
The curve a — g — I — d is an isotherm of pressure P* . When pressure runs from 
a to g, or from g — I the corresponding potential energy changes from U*(n*) 
to U?(n*). The gap U£ - U* g is shown. 

and the interactions between liquid molecules can be described by the same 
interaction potential, such as Sutherland potential. The difference between the 
mean distances of molecules in the two phases results in the difference of molec- 
ular potential energies in the two phases. When a liquid molecule vapors and 
becomes a gas molecule, there is an interaction potential energy gap between 
the molecules in the two phases dU g i. The gap exists only under the SHA and 
the mean distance expansion, in contrast with Mayer's theory in which there 
was no gap. In fact, the partition function (2) was already formally presented 
in some of the previous works, but it is with different meanings and the com- 
parable inferences on phase transition could not be arrived at as no gap existed 
in their discussion. 

The molar latent heat can be expressed simply by 

L = N A (dU gl + Pdv gl ) 

= N A [a{n l -n g ) + P(l/n g -l/n l )}, (13) 

which is equivalent to the Clapeyron's equation[19] and Na is Avogadro's con- 
stant. The meaning of this formula is very clear: the mean distances between 
molecules arc different when the system is in different phases. This difference 
leads to the gap of the intermolecular potential between the two phases, which 
together with the work done due to volume difference between two phases en- 
genders latent heat. Other theories cannot give latent heat also. We define 
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L = L/P c v c and L can be written as 

L = 3W-n;) + P*{±-±). (14) 

which is a function of reduced quantities p*(T*) and n*(T*) only. Fig.lC shows 
the relation between latent heat L and temperature T* . The solid-curve is drawn 
from experimental data of Argon[22]. The dash-curve is the latent heat of VDW 
fluid from Eq.(14). The consistencies between the latent heat calculated by the 
Eq.(13) and that by the Clapeyron's equation provide a strong support to the 
new theory in this letter. 

From Figurel we find that the densities of coexistent vapor-liquid, the vapor 
pressure and the latent heat of VDW fluid are qualitatively correct, but quan- 
titatively incorrect, when compared with the experimental data[20-22]. This 
inconsistency is due to the incorrect potential we have chosen. 

Using the effective potential u(r) — ~eT s (a/r) M (referred to as modified 
VDW fluid) with S=-0.49 and M=2.83 (in Mayer's and Eyring's works, we 
should have M > 3 to keep convergence, but there is no such restriction in 
mean distance expansion), instead of the attractive potential u(r) in Sutherland 
potential, and following the above standard calculation procedure, we obtain the 
better results(see dot-line in Fig.l). 

Now let us review what is new in this work. The starting point of this work 
is the SHA. The difference between the SHA and the previous works, such as 
Mayer's and Eyring's work, is the following: a). The EHSE SHA claimed that 
the number of the particles in a volume d 3 r is dN = N /VA-Kr 2 dr whereas the 
SHA postulates that each particle occupies a volume V/N. b). Mayer and 
Eyring assumed that the distributions of particles are nonexclusive, i.e. the 
distribution of an arbitrary particle is not affected by other particles. On the 
contrary, the SHA postulates that the distributions of particles are exclusive, i.e. 
each volume element V/N can, but can only lodge one particle averagely. This is 
a fundamental difference between mean distance expansion and Mayer's cluster 
expansion, c). The intermolecular interactions in Mayer's and Eyring's calcu- 
lations are independent of the mean distance r, meaning that when the mean 
distance or the density varies, the intermolecular interactions do not change. 
As a result, based on its interactions, we cannot distinguish which phase the 
system is in: vapor or water. But in this work, when the density of the system 
goes from n g to n t , or from n a to n g ,the interaction will change as the mean 
distance varies. 

The new idea in this work is that the total intermolecular potential energy 
and partition function for a classical system in equilibrium are relative to the 
average distance r of molecules (or density). Starting from it, we obtain the 
energy gap between two phases, and a new explanation for the first order gas- 
liquid phase transition is given as follows: first order gas-liquid phase transition 
does not correspond to singularity of thermodynamic functions. In Fig. 2 if a 
system runs from a — g on the isotherm a — g — I — d, the corresponding potential 
energy increases as density changes. At g or I, there are two possible phases wiith 
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different intermolecular potential. These different fluid phases correspond to 
different solutions of Eqs. P (T,n g ) = P (T, ni) and [i(T,n g ) = fi (T, ni)., which 
form two branches of the solutions. For each branch fi t (T,m) or fi g (T,n g ), 
either function /z z (T, m) or fj, g (T, n g ) are analytic (see ref.[19] for detail). When 
T = T c , we have ni = n g = n c and the gap dU g i — > 0. Here there is no request 
for the thermodynamic limit in strict mathematical sense. 
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